Density and Correlation functions of vortex and saddle points in open billiard systems 
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We present microwave measurements for the density and spatial correlation of current critical 
points in an open billiard system, and compare them with the predictions of the Random Wave 
Model (RWM). In particular, due to a novel improvement of the experimental set-up, we determine 
experimentally the spatial correlation of saddle points of the current field. An asymptotic expression 
for the vortex-saddle and saddle-saddle correlation functions based on the RWM is derived, with 
experiment and theory agreeing well. We also derive an expression for the density of saddle points 
in the presence of a straight boundary with general mixed boundary conditions in the RWM, and 
compare with experimental measurements of the vortex and saddle density in the vicinity of a 
straight wall satisfying Dirichlet conditions. 
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I. INTRODUCTION 

The use of Gaussian random functions to describe the 
spatial structure of complex physical systems has had a 
wide range of success, originating with Rice's description 
of the random currents of shot noise [l[ and Longuct- 
Higgins' description of random water waves 0,01 , more 
recently in such diverse fields as sound waves and acous- 
tics H, H[ , turbulence @ , optical speckle patterns^ , and 



the cosmic microwave background fluctuations 

In the realm of quantum wave physics, the same uni- 
versality of the amplitude fluctuations has been conjec- 
tured in the spatial patterns of eigenfunctions in systems 
with classical (ray) chaotic dynamics Q . The analogy be- 
tween the equations of a non-interacting two-dimensional 
electron gas and the electromagnetic modes of a mi- 
crowave cavity 0, [U El EE 03 (see Refs. M, M for re- 
views) allows a unified treatment in the language of quan- 
tum billiards, namely as solutions of the two-dimensional 
Helmholtz equation 



(i) 



with wavenumber k, energy fc 2 , where f = (x,y). We 
therefore may study the properties of electron wavefunc- 
tions, which are difficult to access experimentally, by 
means of measurements using our microwave experimen- 
tal setup. Previously, insights from this analogy have 
had a strong impact on the theoretical study of coherent 
effects on electronic systems in the mesoscopic regime, 
where the spatial correlations of the electronic wavefunc- 
tion are, besides the fluctuations of the energy spectra, 
the source of mesoscopic reproducible fluctuations (see 
Refs. [U M 

for recent examples). 



Concerning fundamental questions, the measurements 
of complicated statistical measures (namely, averages 
over the experimentally constructed eigenfunctions of 
complicated functionals) are very stringent probes for 
the statistical assumptions upon which theoretical mod- 
els of chaotic wavefunctions are based. The primary such 
model is the so-called Random Wave Model (RWM), pro- 
posed by Berry Q, based on the isotropic 2-dimensional 
random waves studied by Longuet-Higgins @ . The basic 
RWM is a statistically stationary isotropic solution of the 
two-dimensional Helmholtz equation (time-independent 
Schrodinger equation), statistically invariant to transla- 
tion and rotation. 

The RWM predicts that the spatial fluctuations of 
cigenfuctions are Gaussian distributed, and this gives rise 
to characteristic morphological features (e.g. Refs. [20l . 
l2l[ ). So far, and largely due to the experimental possibil- 
ity given by the microwave measurements, the assump- 
tion of Gaussian statistics has successfully passed very 
demanding tests. To mention just two examples where 
very complicated functionals of the measured eigenfunc- 
tions are required and the RWM provide excellent results 
we have the distribution of current [22| , the intensity dis- 
tribution in the transition from closed to open billiards 
[23[ and the distribution of quantum stress tensor [24| 
(see Ref . [l6| for a recent review) . 

In this paper we address a different type of functional, 
based on the nodal properties of complex chaotic wave- 
functions, and show again how the assumption of Gaus- 
sian statistics is strongly supported by the experimental 
results. 

The features of the random eigenfunctions we study 
here are the critical points of the current density (here- 
after current) associated with the complex wave tp. The 
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current density is defined by 

j(r) = Im^(f)*VV>(r). (2) 

In quantum- mechanical systems j(r) is representing the 
probability current density at position r. In quasi-two- 
dimensional electromagnetic microwave billiards there is 
a one-to-one correspondence of j(r) to the Poynting vec- 
tor [23. 

Since ip is assumed to satisfy the two-dimensional 
Helmholtz equation, the points where j — are of two 
types: vortices of the flow (also known as circulations, 
wave dislocations, nodal points and phase singularities 
[U, SI H, [H, [27[ ), where V = and about which 
the current swirls in a counterclockwise (+) or clock- 
wise (— ) sense, and saddle points (stagnation points), 
which are also saddle points of the phase arg-0, and hy- 
perbolic points in the current flow. (The existence of 
phase extrema - maxima or minima - is prohibited by 
the Helmholtz equation [2lJ •) The topological Poincare 
index of these types of points, describing the number of 
turns of j in a small circuit of the critical point, is +1 
for vortices (regardless of the sense of circulation), and 
— 1 for saddle points. General arguments based on sta- 
tistical isotropy demonstrate that there can be no net 
accumulation of topological charge, cither in the sign of 
the vortices, or in the Poincare index. Therefore, the 
densities of positive and negative vortices must be equal, 
and the bulk vortex density must equal the bulk saddle 
density; calculations based on the RWM [H, HE give 
this vortex density as k 2 /Air. Knowledge of the positions 
of the critical points of the current vector field provides 
a skeleton on which the rest of the flow field is based. 

An important deviation to the bulk RWM is caused 
by the presence of boundaries. The interplay between 
spatial confinement and Gaussian fluctuations is by no 
means trivial, and it has even been claimed that in con- 
fined systems Gaussian fluctuations are valid only over 
very short distances [28|. Substantial progress has been 
made recently in modifying the basic RWM to include 
boundary effects [H, [iEjp, [H, H3, H3, HI] , supported 
by numerical evidence |3(|. and there appears to be no 
reason to believe that the Gaussian assumption fails in 
the vicinity of a system's boundary. Here, we present 
to our knowledge the first experimental measurements 
demonstrating the validity of boundary-adapted RWMs, 
based on Dirichlet conditions on an infinite straight wall. 

The intensity distributions in a two-dimensional 
chaotic microwave cavity, as shown in Fig. [1] are well- 
understood [23[ , although the distribution of current crit- 
ical points has not previously been studied in detail. 
In addition to the density fluctuations of vortices and 
saddle points against a Dirichlet boundary, we measure 
the vortex- vortex correlation function (including the case 
signed by circulation), and the vortex-saddle and saddle- 
saddle correlation functions, comparing against predic- 
tions of the RWM. 

Although some of the theoretical predictions we com- 
pare with have been derived before (such as the vortex- 




FIG. 1: (color online) Morphology of typical complex wave- 
function tp in the open quantum billiard of our experiments, 
(a) Modulus (intensity) \tp\ 2 ', (b) Current flow Imi/>*Vi/>; (c) 
Blowup of (b), demonstrating the critical points we study. 
The vortices are marked by dots, colored according to sense: 
counterclockwise (red), clockwise (green). The saddle points 
are marked by crosses. Our points of measurement are the 
crossing points in the background grid. 

vortex correlation functions 0, [2(| and vortex density 
fluctuations against a straight wall [2^, H(|), others are 
new. In particular, we derive the density fluctuations 
of saddle points in the presence of an infinite straight 
boundary on which the wave satisfies mixed (Robin) 
boundary conditions, and, although we were unable to 
derive exact saddle correlation functions analytically, we 
have found large-r asymptotic approximations to these 
functions. 

The paper is organized as follows. In section |TT] we 
explain the experimental set-up and techniques used to 
locate the critical points in the microwave cavity. General 
definitions and properties of critical points of the current 
associated with the Helmholtz equation, as well as the 
expressions for densities and correlations, are described 
in section ITTIl This is followed by section ITVl which de- 
scribes our RWM calculations: the model is introduced 
in subsection IIV Al and details follow for bulk correlation 
functions (subsection IIV Bp and densities near a straight 
boundary satisfying Robin conditions ( subsection IIV C[) . 
The experimental results are compared with the theoret- 
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ical predictions in section fVl 



II. EXPERIMENTAL SETUP 

We here report on the measurement of correlations 
functions of saddle and vortex points in an open billiard 
system including effects of the boundaries. The basic 
principles of the experiment can be found in Ref. l37l We 
used a rounded rectangular cavity (21 cmx 16 cm) cou- 
pled to two wave guides of width 3 cm with a cut-off 
frequency at ut = 5 GHz. To break the symmetry and 
to block direct transport, two triangular obstacles with 
a length of 12 cm and a height of 1 cm were placed in the 
resonator. Absorbers were placed at the end of the leads 
to avoid reflection. We scanned the billiard on a square 
grid of 2.5 mm with a movable antenna A\ and measured 
transmission S\2 in the range of 4- 18 GHz from a fixed 
antenna A2 in the end of the right lead. The fixed an- 
tenna had a metallic core of diameter 1 mm and a Teflon 
coating while the probe antenna A\ was a thin wire of 
diameter 0.2 mm to minimize the leakage current. The 
lengths of the antenna A\ and antenna A2 were 4 and 
5 mm respectively. 

For microwave frequencies v < c/2d = 18.75 GHz, 
where c is the velocity of light and d is the resonator 
height, the billiard is quasi- two-dimensional. In this 
regime there is an exact correspondence between elec- 
trodynamics and quantum mechanics, where the compo- 
nent of the electric field perpendicular to the plane of 
the microwave billiard E z corresponds to the quantum- 
mechanical wave function ip. Figure [1] shows typical in- 
tensity and current patterns thus obtained. Additionally 
a zoom of the flow pattern is shown to visualize the struc- 
ture and showing the different types of critical points in 
the flow. 

We previously have reported results on vortex pair cor- 
relation functions and nearest neighbor distance distribu- 
tions [38l . [HsJ for vortex points. Through improvements 
in data analysis, we have been able to study the sad- 
dle points of the current in addition to the vortices, and 
significantly reduce the effects of noise in our measured 
correlation functions. These improvements are sketched 
in the following. We have increased the spatial resolu- 
tion by a factor of two compared to our previous mea- 
surements before and we now use additionally a bilinear 
interpolation method for the individual components of 
the flow. Using the bilinear interpolation we estimate 
the nodal lines of the individual flow components, en- 
abling us to get the exact position for vortex and saddle 
points within this approximation. The extraction of crit- 
ical point locations is now fully automatic, allowing large 
samples of data to be analyzed. One can see the effect 
of these improvements in the better results for the pair 
correlation functions (good agreements also for small kr) 
and the charge correlation function which we can present 
here for the first time. 



III. CURRENT CRITICAL POINT DENSITIES 

In this section and following, we assume that the com- 
plex wavefunction ip(r) = £(f) + vrj(f) has no particular 
symmetries or properties, beyond satisfying Eq. |T]). The 
current j, from Eq. ([2]), can therefore be written 

3 = (€Vx - Vtx , fry -Vtv)- (3) 

j = at vortices, where ip = £ = 77 = 0, and at saddle 
points, where £/t) = £ X /Wx = tiy/Vy (if V = 0j there is 
equality between the reciprocals of these terms). 

The quantity which distinguishes vortices from saddles 
is the Jacobian 

J = d X j X 8yjy - dyj X d X jy, (4) 

which is positive at vortices, and negative at saddles. 
Since jl) satisfies Eq. |T]), J separates into two contribu- 
tions |27j . 

J = Jv- Js, (5) 

where 

Jv = (ZxVv ~ £yVx) 2 , (6) 

HtV*v - V^y) 2 - (7) 

Obviously, J v = at saddle points, and J s = at vor- 
tices. This fact, combined with positivc-definiteness of 
the two parts of the Jacobian, implies that \J\ = J y at 
vortices, and \ J\ = J s at saddle points. 

These quantities can be used to define functions which 
find critical points (vortices or saddles) at position f. The 
density of critical points, with a unit (5-function at each 
zero point of j , is given by 

£ crit (f)EE J2 s2 ^) = S 2 (j(r))\J'(r)\. (8) 

{r: J(r)=0} 

By the separation of J above, this gives the saddle den- 
sity m, 

D s (r) S 2 (j(rj)J s (r). (9) 
and the vortex density [25|, [2(| H?) , 

D v (f) EE S 2 (j(f))J v (r) = d(0S(v)\$xVy ~ £yTh\- (10) 

The vortex sign (sense of circulation) is given by 

S EE sign(^% - tyflx), (11) 

so removing modulus signs gives the signed vortex den- 
sity. 

The number and location of critical points for a given 
field must be founded by explicitly solving the set of equa- 
tions j(f) = 0. This of course requires the precise knowl- 
edge of the spatial dependence of the particular solution 
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ip(r) in which we are interested. The task of solving 
the Helmholtz equation in cases where the geometry of 
the confinement (transversal section of the waveguide) 
is such that Eq. ([T|) is not separable is usually very de- 
manding. This makes the function-by-function study of 
current morphology almost impossible. 

A suitable way to overcome this complication is to use 
a statistical approach. This idea is based on the strong 
uniformity of the solutions of the Helmholtz equation 
with non-integrable geometries (see for example Fig.[Tj), 
indicating that their main properties actually depend on 
far fewer parameters than the full solution itself. We 
therefore consider, instead of a given set of solutions of 
the Helmholtz equation, an ensemble of fields. This en- 
semble will be constructed in such a way that the most 
general and basic properties of the exact solutions are 
respected, in the hope that these general properties suf- 
fice to fix the morphology fluctuations. The ensemble we 
choose is the usual random wave model (RWM) discussed 
in the Introduction, or the boundary-adapted model of 
Refs. [H ML 

Delaying on the appropriate definition of the RWM 
until the next section, we merely write the average over 
the ensemble (■•■)• I n this paper, we compare theoretical 
RWM predictions and experimental measurements of the 
average density fluctuations 

P a (f) = ^{Da(?)), (12) 



k 2 



and 2-point correlations 

9ap(f2,fl) 



4tt 



(D^D^n)), (13) 



where a, are v,s. These expressions have been normal- 
ized against the bulk average vortex density (and saddle 
density) k 2 /4tt (2f| HH, [27} ■ Statistical symmetries in the 
RWMs will mean that the densities p and correlations g 
have simpler functional dependence. 2-point correlation 
functions can also be considered which take topological 
signs into account, such as the vortex topological charge 
correlation function 

(-D v (f 1 ) ( S(ri))i) v (f 2 ) < S(r 2 ))) 



9Q{n,r 2 ) 



(k 2 /4ir) 



(14) 



where the effect of the signum S functions of Eq. (Till) 
is to negate the modulus signs in D v ; the vortices are 
signed by their sense of circulation. Current critical point 
correlation functions can be written down in terms of 
the correlation functions g a p. If the 2-point critical point 
correlation function is denoted g C rit, and gi the function 
signed by Poincare index (positive for vortices, negative 
for saddles), we have 



t(ri,r 2 ) 



gi(n,f 2 ) 



•(n,f 2 ) +g ss (ri,r 2 ) + 2g vs (fx,f 2 )) , 

(15) 



4 {9w(ri,f 2 ) + s ss (ri,f 2 ) - 2flv s (n,r 2 )) ■ 

(16) 



IV. AVERAGE DENSITIES AND 
CORRELATIONS OF CRITICAL POINTS 
WITHIN THE RANDOM WAVE MODEL 

A. RWM: basic definition and field correlations 

The random wave model assumes that the wave field, 
satisfying the Helmholtz equation (|TJ) , is a superposition 
of infinitely many complex plane waves with wavenumbcr 
k with uniformly random directions and phases; the real 
and imaginary parts of the field are therefore assumed to 
be statistically independent. These assumptions, while 
appropriate to our experimental open billiards, which 
lack time reversal symmetry, do not apply in the tran- 
sition from open to closed systems |16fl . The RWM is 
crgodic in the bulk - spatial averages are equivalent to en- 
semble averages. We will describe the boundary-adapted 
random wave model at the end of this subsection. 

The central limit theorem ensures that, in the limit of 
infinitely many superposed plane waves, the probability 
density function of the value of the wave at each point has 
a complex circular Gaussian distribution @,[25|. Further- 
more, the distribution of all derivatives of the field are 
also Gaussian random variables, which may have non- 
vanishing correlations with each other and the original 
field. 

The assumption that the field and its derivatives pos- 
sess multivariate Gaussian statistics implies that, for a 
functional jF[fi], depending on the field and its deriva- 
tives at possibly different points, we have 



1 



v /(2 7 r)"detM 



(17) 



where u is an n-dimcnsional vector consisting of the rel- 
evant Gaussian random fields £(r\), d y rj(f 2 ), etc. appear- 
ing in and M is the n x n matrix of correlations with 
entries Mi ,-, 



Mij = (muj). 



(18) 



Calculating densities of morphological features in the 
RWM is therefore reduced to a Gaussian integral, whose 
difficulty depends on the complexity of the functional T . 
For instance, the average density of vortices or saddles in 
the bulk isotropic random waves can be calculated with 
T = D v or D s from Eqs. @, (JTDJ) ; as discussed previously 
these are known to be equal constants, with value k 2 /Air 
[27l | . In this paper, we concentrate on two specific types 
of functionals T . 

In subsection II V Bl we consider two-point correlations 
g a fj in the bulk isotropic random wave model, where T is 
given by D a (ri)Dp(r 2 ), with a,/3 = v,s. These expres- 
sions are the average densities of vortices or saddles at 
two points, depending only on the scaled distance 



R = k\r 2 - f\\, 



(19) 



by isotropy. These critical point correlation functions de- 
pend only on the 2-point field correlation function, given 
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FIG. 2: (color online) The normalized mean field intensity 
B(Y\ a), plotted against Y for various choices of Robin param- 
eter a: a — 0, i.e. Dirichlet conditions (solid, blue); a = n/4 
(dotted, purple); a = n/2, i.e. Neumann conditions (dashed- 
dotted, yellow); a = — n/4 (dashed, green). The black line is 
at 1 (the asymptotic limit for Y — > oo) 



by 



c{R) = -MWfa)} = = (v(nM?2)). 

(20) 

All 2-point correlation functions of derivatives of the field 
can be expressed in terms of derivatives of C(R). It is 
well known Q that the field correlation function of the 
2-dimcnsional isotropic RWM is given by the Bessel func- 
tion 



C(R) = J (R). 



(21) 



Much of our argument will be based on asymptotic ap- 
proximations for large R, in which, to leading order, 



C{R) 



■ cos(i? - 7r/4). 



(22) 



The second type of functional we consider is based on 
the vortex and saddle densities D V ,D S in the so-called 
boundary adapted random wave model. In this model, 
the wave with y > is assumed to satisfy a homogeneous 
boundary condition along the infinite straight line y = 
0. As with the correlation function above, we will use 
dimensionless Cartesian coordinates 



X = kx, Y = ky. 



(23) 



Although our experimental data are for Dirichlet condi- 
tions only ("0(A, 0) = 0), our discussion will be framed in 
terms of the most general boundary conditions, namely 
mixed (Robin) conditions 



ip(X,0) cosa + dYip(X,0)sina = 0, 



(24) 



where a is a constant with < a < n. Dirichlet condi- 
tions correspond to a = 0, and Neumann to a = tt/2. 



Berry and Ishio [30( constructed a natural RWM sat- 
isfying Eq. (j2~4"|) , and considered the vortex density as 
a function of distance Y from the boundary, and gen- 
eral a, generalizing previous work [29j | for Dirichlet and 
Neumann conditions. We will calculate the correspond- 
ing y-dependent saddle point density in subscction llV CI 
The critical point density calculations depend on the In- 
dependent, 1-point quadratic field correlation 



(25) 



For general a, the function B{Y\ a) for the mixed bound- 
ary condition RWM of Ref. [30 cannot be expressed in 
terms of elementary functions, although it has a straight- 
forward integral representation: 

B(r;o) = l-| £ 2 d9 (26) 

(1 - t sin 2 6) cos(2y sin 6) + 2t sin 9 sin(2Y sin 6) 
1 + t sin 2 9 

with r = tana. For Dirichlet and Neumann conditions, 
B can be expressed in terms of the Bessel function Jo, 



(27) 



It is straightforward to find an asymptotic approximation 
for B(Y;a), 



B(Y; a) Y ^ 1 



'■kY 



cos(2y - 2a- tt/4), (28) 



consistent with Eq. (|2"7| . The field intensity fluctuation 
B(Y; a) is plotted as a function of Y for various choices 
of a in Fig. [2] We mention that the result in Eq. ([28]) 
is also found by using the semiclassical approximation 
for the two-point correlation function in the presence of 
boundaries with Robin boundary conditions. In this case, 
the parameter a enters through the semiclassical phase, 
as explained in [35| . 



B. Spatial correlations of current vortices and 
saddles: bulk results 

Correlations of vortices and related objects have been 
the subject of much study in the isotropic RWM. In 
particular, the signed vortex-vortex correlation function 
9q{R) defined in Eq. 11141). i s known to have a remarkably 
simple form [H, Sfl EJll [H, 0, H| 



' 9Q(i?) = ldl? 



darcsin(J (i?)) 



n 2 



dR 



(29) 



This equation (with Jo replaced by a suitable 2-point 
function) holds for general isotropic Gaussian random 
fields, not just the Bessel-correlatcd random wave model. 
On account of the isotropy of the distribution in the 
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phase of the field ip, gQ satisfies the 'topological charge 
screening relation' (i gno r ing the self-interaction at R = 

o) m, ra, sn, m, us, S3, eh 

±J dRRg Q (R) = -l. (30) 

For R > 1, g Q (R) ~ 8cos(2R)/ttR 2 . The oscillation pe- 
riod of gQ^R) is twice that of the correlation function 
C(R) since there are two nodes per oscillation of J . 

However, the unsigned correlations we emphasize here, 
including saddle correlations, do not have such a simple 
form. The vortex- vortex correlation function g vv (R) was 
computed exactly in Refs. l25l . I2H and is written down 
in Ref. \M Eqs. (32)- (35) (also see Ref. S3 p. 83) as a 
complicated function involving various elliptic integrals; 
using similar methods involving computer algebra, our 
attempts to extract g vs (R) and g ss (R) analytically were 
unsuccessful. We therefore developed a strategy, based 
on an asymptotic expansion of the correlation matrix 
through the asymptotic expansion of the Besscl function, 
as indicated in the subsection llV Al to extract asymptotic 
approximations for these correlation functions (including 
9w(R)) to compare with the experimental data. Details 
are given in appendix [X] Our asymptotic argument is 
similar to the asymptotic 2-point correlation function for 
gradient saddles in computed in the real RWM [48[ , based 
on functional differentiation. 

In order to compute the saddle-saddle spatial correla- 
tion we need to consider as degrees of freedom the field 
and its derivatives at two different points, fi, r<i. The ba- 
sic idea of the asymptotic method is that in this case the 
correlation matrix in Eq. (|17[) can be unambiguously sep- 
arated into two contributions. On one side, we have the 
correlations of the degrees of freedom at the same point, 
which is of course independent of R = k\f\ — fa|. All the 
spatial dependence of the correlation is encoded in the 
elements of the correlation matrix relating the fields at 
different points. 

Once this separation of the spatial dependence of the 
correlation matrix is realized, the asymptotic expansion 
of the two-point correlation function Eq. (|2"2")) naturally 
leads to an expansion of any arbitrary average in powers 
of the small parameter 1 / \/R. Obviously, this expansion 
will fail for short distances (comparison with the experi- 
mental results show that in practice the approximation is 
good already for R > 3). Omitting further details of the 
derivation, the relevant asymptotic approximations to or- 
der 0(R~ 1 ) of the RWM vortex- vortex, vortex-saddle, 
and saddle-saddle pair correlations are 



g vv (R) " 




4 sin 2R 


(31) 


- H 


nR ' 


9vs(R) - 




4sin2i? 


(32) 


- 1 - 


nR ' 


9ss(R) " 




4sin2i? 


(33) 


- H 


nR 



Fig . is a plot of <? vv (R) , computed analytically 
[46j , and from the asymptotic form (p?T]) . The 2- 



1 Cj I — i — r-p — | — i — i — i — | — i — i — i — | — i — i — i — | — i — i — i — | — i — i — i — | — i — i — i — | — r 




— 0.5 I i i i I i i i I i i i I i i i I i i i I i i i I i i i I l 

2 4 6 8 10 12 14 
R 



FIG. 3: (color online) Analytic 2-point correlation functions 
g vv (R) (red) and Qq{R) (green), plotted against R. The black 
line is the asymptotic value of 1 for g vv (R) as R — > oo. The 
analytic functions compare well with their asymptotic approx- 
imations (dashed curves) for R > 4. 



point correlation functions for vortices of the same sign 
9++(R) = .9 — (R) = |(sVv(-R) + 9q{R)) and opposite 
sign g + -(R) = ^(g vv (R) — gg(R)) oscillate in phase plj | 
since gQ ~ 0{R~ 2 ) decays more swiftly than g vv (R). 
These equations demonstrate that critical points, that is 
vortices (whose Poincare index is +1) and saddles (with 
index —1) oscillate out of phase. 

Eqs. (|31|) - ((33|) can be used to estimate asymptoti- 
cally the critical point 2-point functions g C mt(R), and its 
Poincare index-signed analogue gi (R) : 

gcrit(R) ~ i + 0(R- 2 ), (34) 
9l (R) ~ 4sin(2i?)/7ri?. (35) 

gi(R) decays rather slowly, in contrast to the long range 
correlations of topological charges of other RWMs, such 
as critical points of the gradient in the real RWM [Z3L l44l| . 
which decay to leading order like 0(R~ 3 ). It would there- 
fore be interesting to establish whether Poincare index 
satisfies a screening relation analogous to Eq. (|3"0|) . since 
convergence in the integral is marginal. 



C. Densities of critical points near mixed 
boundaries 

In the boundary-adapted RWM, the density of vortices 
and saddles oscillates with distance Y from the boundary, 
just as the square field does (e.g. Eq. (|28|) ). However, 
with the Robin RWM of Ref. [30, the form of the vortex 
density function depends only on the function B(Y;a) 
of Eq. (f2"5"]): the result of the actual Gaussian integral is 
independent of the value of parameter a. The entries of 
the correlation matrix, of course, do depend on a, and all 
may be written as linear combinations of B(Y; a) and its 
derivatives with respect to Y. 
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Therefore, the mean vortex density in the boundary- 
adapted RWM, is [H[ 



y/±B + B»-\y/B{2 + B»)-B'* 
Pv(Y; a) = ^j- 2 , (36) 

where dependence of B on Y and a is suppressed and the 
prime denotes the partial derivative with respect to Y. 
The density is normalized with respect to the bulk den- 
sity, so limy^oo p v (Y; a) = 1. This equation is the same 
as Ref. HiJ Eq. (40), with correlation matrix elements 
replaced by appropriate functions of B. 

Since the saddle density D s (r) of Eq. © does not in- 
volve modulus signs, the calculation of the average saddle 
density in the boundary-adapted RWM uses straightfor- 
ward Gaussian integration techniques, as outlined in Ap- 
pendix [B] The resulting density is 



Ps(Y;a) 



B 1 ' 2 



(45 + B" - 2) 3 / 2 (B(2 + B") - B' 2 ) 3 / 2 
x (16 - 645 + 645 2 + 165' 2 - MBB' 2 
+ 165' 4 - 165" + 645 2 5" - 3255' 2 5" 
+ 16BB" 2 + 165 2 5" 2 - 45' 2 5" 2 
+45" 3 - B" 4 + SB'B'" - 3255'5"' 
+8B' 3 B"' - 8B'B"B"' + 2B'B" 2 B"' 



-AB 2 B"' 2 



B' 2 B" 12 - 2BB"B"' 2 



-855"" 
-8BB' 2 B"" 



-2B' 2 B"B"" 



165 2 5"" 



45' 2 5"" 



85 2 5"5"" 
2BB" 2 B"") 



(37) 



Asymptotically, for 7> 1 we find 



p v (Y;a) ~ 1 + 



2cos(2(y-o)-f) l + 5sin(4(y-o)) 



4ttY 




0.300 



0.310 



0.330 



0.340 



(38) 



Ps(Y 



2cos(2(y-a)-5) l + 5sin(4(y-o)) 
; a ) ~ 1 1 T^Cr 



47ry 



(39) 



(Eq. (|55]l was demonstrated in Ref. [30|). Therefore the 
leading order oscillations in vortex and saddle densities 
are exactly out of phase, as in the 2-point correlation 
functions discussed in the previous subsection. Also, 
as with the correlation functions, the oscillations have 
twice the periodicity of the underlying correlation func- 
tion (again, as nodes occur with double the frequency of 
a sinusoidal wave). The mean saddle density for several 
choices of a is plotted in Fig. 2^. 

As shown in Ref. H3, when a <C 1, there is a large 
additional peak in the vortex density for small R. Robin 
boundary conditions in this regime are known to have un- 
usual properties, such as admitting negative energy solu- 
tions [49l | , and diverging momenta on the boundary [50l | . 
When a is small, the peak occurs in the neighborhood 



FIG. 4: (color online) Density oscillations of critical points as 
a function of distance Y from a wall satisfying mixed Robin 
conditions: (a) p v (Y;a); (b) p s (Y;a). The colors represent 
the same choices of a as in Fig. [2] the black line is at 1 (the 
asymptotic limit for Y — » oo); the two densities are clearly out 
of phase for Y > 4. (c) Peak in the densities for a = 7r/10 : 
p v (red curve), p s (green curve), with the dashed lines the 
small- a forms of Eqs. (gDJ, (gTJ. 



form [3C 



p v (Y = a + a 3 e; a) 



12^1 + 4(6e-l) 2 
a 3 (l + (12e- l) 2 ) 3 / 2 



(40) 



with a peak near Y = a + a 3 /12 of approximately 
24^/a 3 . Analysis of Eq. {ST} reveals that the saddles, 
too, have a peak for small a, and an analogous argument 
as for vortices gives for a <C 1, 



y 



ea , where p v (Y;a) has the skewed-Lorcntzian 



p v (Y = a + a 3 e; a) 



12^2 + 24e(6e - 1) 
a 3 (l + 4(6e-l) 2 ) 3 / 2 



(41) 
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FIG. 5: (color online) Vortex pair correlation g vv (R) and 
charge correlation function gci(R), in (a) the low frequency 
regime (5 GHz < v < 9 GHz), and (b) the higher frequency 
regime (15 GHz < v < 18.6 GHz). The solid lines corre- 
spond to the analytic prediction of RWM discussed in sub- 
section [TV£] using the exact formula for g vv (R) rather than 
the asymptotic form. 



that is, a peak of almost the same shape and magnitude 
as for vortices, but with maximum near Y = a + a 3 /6. 
This peak ensures that small a does not give rise to an 
accumulation of total Poincare index near the boundary. 
Plots of the p v (Y;a),p s (Y;a) for a = n/10 are shown 
in[4j;, along with the corresponding approximations. 



V. COMPARISON BETWEEN EXPERIMENT 
AND RWM PREDICTIONS 

In this section we compare the results of our microwave 
billiard experiment, outlined in section [Til with the theo- 
retical predictions of the isotropic and boundary-adapted 
RWMs described in Section HVl Since distances R,Y are 
scaled with respect to k as in Eqs. §T§\i . (|23[) . the ex- 
perimental results for different frequencies v have been 
superimposed, improving the statistics. 

Fig. [3] shows the comparison between experiment and 
theory for the vortex-vortex pair correlation functions, 
both signed (Eq. and unsigned (Ref. 51 Eqs. (32)- 
(35)), in two different frequency regimes. Because of the 
limited resolution due to the measurement grid in the 




n n h i i i - 

10 20 30 40 

R 




10 20 30 40 

R 



FIG. 6: (color online) Pair correlations involving saddle 
points: (a) saddle-saddle correlation function g ss (R); (b) 
vortex-saddle function g vs (R). Experimental data is plotted 
against the asymptotic forms of Eqs. (|32[) . (|33ll . 



low frequency regime the correlation function could only 
be determined reliably for small kr, see Fig. [5jt. But in 
this regime the system size is comparable to the vortex 
spacing, leading to an influence of the boundary on the 
measured correlation function and restricting the observ- 
able R range up to i? max = kL ~ 15, where L corresponds 
to a typical length of the system. Due to the boundary 
effects the oscillation period of the experimental results 
is shorter than the theoretical one for the bulk statistics 
[Sll l53 j . The results for the higher frequency regime, on 
the other hand, shown in Fig. [5Jd agrees perfectly with 
theory for large R, but fails for small R because of the 
mentioned limited resolution in the measurement. Ex- 
perimental results for the vortex pair correlation function 
have been published already previously, though with a by 
far poorer statistics (39j . All other quantities shown in 
this section have not been published previously. 

Results for the saddle pair correlation function g ss (R) 
and the vortex-saddle function g vs (R) are plotted in 
Fig. [BJi, with experimental data plotted against the 
asymptotic formulas (|3"3")h (f3"2"]) . As with the vortex cor- 
relation functions, the agreement between theory and ex- 
periment is very good for a wide range of R; since the 
theoretical formulas arc asymptotic, we do not have a 
theory to compare with the data for small R. Fig. [SJd 
shows the total and signed critical point correlation func- 
tions g cx it(R) and gi(R), again with good agreement. 

Experimental measurements of the average vortex and 
saddle density fluctuations against a straight boundary 
satisfying Dirichlct conditions are shown in Fig. vortex 
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FIG. 7: (color online) Critical point density fluctuations as a 
function of scaled distance Y from a straight wall satisfying 
Dirichlet boundary conditions: (a) vortex density; (b) saddle 
density. Experimental data is plotted against the analytic 
forms of Eqs. (|36p . (|37|l . for Dirichlet boundary conditions. 



fluctuations in the upper panel (with theoretical density 
from Eq. ([55]) ). the lower panel the saddle fluctuations 
(with theoretical plot from Eq. (|37[) ). The agreement 
between theory and experiment is excellent. 

In summary, by applying a higher grid resolution and a 
bilinear interpolation technique the flow pattern through 
an open microwave billiard could be resolved by nearly 
one order of magnitude better as in previous experiments. 
This allowed the determination of various distribution 
and correlation functions for the critical points in the 
flows, vortices and saddles, which had been inaccessible 
hitherto. 

From the theoretical side, we have presented the uni- 
versal predictions of the Random Wave Model for the pair 
correlations of the current's critical points in the bulk, 
showing excellent agreement with the measurements. Al- 
though some of the results can be obtained in closed 
form, an asymptotic method valid for large separations is 
necessary to construct some important correlations. By 
power counting of the characteristic decay of correlations 
with respect to the separation screening can be explicitly 
checked, leading to a surprising behavior of the Poincare 
index, which is still to be explained but it is fully sup- 
ported by the experimental results. 



We also contrast for the first time the basic assump- 
tions of the boundary-adapted Random Wave Model 
against experimental results. The influence of the bound- 
ary showed up first in an oscillatory behavior in the den- 
sity of vortices and saddles close to the wall. 



APPENDIX A: ASYMPTOTIC PAIR 
CORRELATIONS 

In this section we sketch the method we use to calculate 
the asymptotic (large R) approximation to the RWM av- 
erage of an arbitrary functional T depending on the field 
and its derivatives at points fi,f*2 with k\f\ — r*2| = R 
the scaled distance. 

We start with the exact expression for the Gaussian 
average 



1 



x /(27r)"+ m detM J- 
where the vector 



^"[u]e"5 



S-M 



(Al) 



(A2) 



comprises all the relevant degrees of freedom ~ 
(u± , . . . , ui 1 ') at position r\ and = {u{ \ . . . , u$) 
at position r*2. With this definition, the correlation ma- 
trix has a natural block form 



M = 



M (2,l) M (2,2) 



(A3) 



where [M^)] = (uf ] uf ] ). The key step of the 
method is to observe that the only dependence of the 
average (J 7 ) on the scaled distance R comes from the off- 
diagonal blocks, and the known asymptotic expansion of 
M^' 2 )(i?) will imply an asymptotic expansion of (J 7 ). 
In order to follow this program, we found convenient to 
switch to the Fourier representation of the probability 
distribution in Eq. (|A1[) 



iiT-M 1 -u 



^(27r) Tl+Tn detM 



d «^(i) d ™^(2) e i(» (1) ^ (2) )-(« (1) ,s (2) ) 

(A4) 



x c _I( t z;< 1 )^( 2 )).M-(«7( 1 ',tu< 2) ) 



Due to the asymptotic form of the off-diagonal block 
M^ 1,2 ) ~ 1/VR, the last line in Eq. ([A"3]) can be writ- 
ten as 



g-i^ 1 ',^ 2 ')-]^-^ 1 ',^ 2 ') _ 

^i^D.Mf 1 ' 1 '.^ 1 '^!^ 2 '-]^' 2 ' 2 '.^ 2 ' x 



(A5) 



-O ( l/i? 3/2 



We note that, given the particular form of the asymptotic 
expansion of the Bessel function, the very first term of 
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the asymptotic expansion of the off-diagonal blocks not 
only gives the leading order term in 1/y/R of the average, 
but also the subleading one of order 1/R. Beyond the 
subleading order, higher order terms of the average re- 
quire higher order terms in the expansion of M^ 1 ' 2 ) . This 
has to be consider a very special property of the Bessel- 
correlated RWM with its characteristic slow (~ 1/y/R) 
decay of correlations. 

Substitution of Eq. (fA"5j) into Eq. ifAl]) results in 
the asymptotic expansion of the probability distribution 
which in turn leads to the sought asymptotic expansion 
of the average in Eq. (|A1[) . The calculations are simple 
but tedious, so we merely quote the result. Denoting 



A = 



M (l,2) 



M 



(2) 



and introducing the tensors 



^a,P) _ /„(<*) T \ 



(A6) 



(A7) 



where (. ..)o indicates the average in Eq. (|AT]) with 
M (i,2) = M (2,i) = 0) we get finally 



{F) = Fo + TrAF 



(2,1) 



-F TrAM (2 - 2) A T M (1 ' 1 
2 

1 



(A8) 



■RF^AM^A' + TrF^AHV^^A 



(2,2) 



A r F^ ia) +0 (l/R 3/2 



+ 2 TrAF o 

Our result Eq. (|A8|) allow us to calculate the leading 
(~ 1/y/R) and subleading (~ 1/R) contributions to any 
correlation in terms of the M^ 1 ' 2 ) = M^ 2 ' 1 ) = (uncor- 
rected) results. It also provides the large R asymptotics 
to one-point functionals at points separated a distance 
2R from an infinite straight boundary. 



APPENDIX B: DERIVATION OF THE 
T-DEPENDENT SADDLE DENSITY EQ. ([37)) 

The saddle density at scaled distance Y from a straight 
boundary can be calculated by standard methods of 
Gaussian integration (as used, for instance, in Refs. l2ll . 
EH), although the details are rather tedious and only 
outlined here. Normalized by the bulk density, the In- 
dependent saddle density is 



P S (Y) 



MDy(Y)) = ~ J d 2 f(cxp(ij ■ t)J s ). (Bl) 



Since J7s is a sum of terms, linearity of the average sim- 
plifies the Gaussian average to calculations of the form 

Ax/3 = ((fiVaP - V£ a p) 2 cxp(ij ■ *)}, (B2) 
where a(5 = XX, YY or XY. 



Each A a p is an average over an 8-dimcnsional Gaussian 
random vector 

W= tap, V>VX,TlY,Tl a i3}. (B3) 

The terms in £ and rj are uncorrclatcd. Denoting cither 
by v, and suppressing y-dependence, the relevant non- 



vanishing correlations follow from Eq. (f2l>|). B 


(vv Y y) = 


1 1 nil 

2 ~~ i B ' 


(4) = 


-\+B + \B", 


(4) = 


\ + \B", 


(wy) = 


2-B', 


(v Y V X x) = 


in' i n'" 
--B - -gB , 


(v y v Y y) = 


1 rtlll 

8 B > 


(v 2 xx) = 


-I + B + 1B" + ±B"" 


(v 2 yy) = 


3 , 1 ryim 

s + TE B ' 


(v 2 xy) = 


1 i 1 Tyii i 1 nim 

& + i B +Te B > 



(B4) 



where the prime denotes the partial derivative with re- 
spect to Y. We denote the appropriate correlation matrix 
for u by Mg. 

The scalar product in the exponent in (|B2[) can be 
written as a symmetric quadratic form in u, 



j ■ t = hu ■ T • u, 



(B5) 



where T depends on t\ , ti . 

Therefore, defining the matrix S = (M^ 1 +iT) _1 , and 
Q = (£, r lap — vtap ) , it is straightforward to see 

^ = (2,)V 1 detM, / d8 ^ 2eXP( -^- S " 1 -^ 



' detH 
dctM 



-[A 2 exp(-±w-3.w)], =0 (B6) 



where w is a Fourier dual to u, and A is the quadratic 
form of Fourier derivatives corresponding to Q. The final 
step in Gaussian integration by parts reduces to a mul- 
tilinear combination of entries of S. Each A a p can now 
be integrated with respect to t. In terms of the original 
correlations, the final result is 



Ps(Y) 



(4)3/ 2(jB(4) _ (^5)375 [( B ^> - <^> 2 ) 

«4>«4*> + «y) + 2(4y» - 2(vyvxx) 2 ) 
-(v 2 x)(vy)((vx) 2 + (wyy) 2 ) 
-B(v 2 x )((v y vxx) 2 + (vyvyy) 2 ) (B7) 
-2(v%)(vvy)((v x ){vyVxx) + (w Y y)(vyvyy))] ■ 



Eq. (|37p follows from this expression with the appropriate 
substitutions from Eqs. (|B4[) . 
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